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Achieving a precise description of multi-parton final states is crucial for many analyses at LHC. In this contri- 
bution we review the main features of the HELAC-NLD system for NLO QCD calculations. As a case study, NLO 
QCD corrections for ti + 2 jet production at LHC are illustrated and discussed. 



1. Introduction 

With the successful start of collisions at 7 TeV 
a few months ago, the Large Hadron Collider 
has put another milestone towards a thorough 
exploration of the Terascale. The large energy 
available will open many multi-particle channels 
that have being widely investigated pQ. In view 
of a correct interpretation of the signals of new 
physics which might be extracted from data, it 
is of considerable interest to reduce our theoret- 
ical uncertainty for the physical processes under 
study, especially when large QCD backgrounds 
are involved. In this respect, the need of Next- 
to-Leading- Order (NLO) corrections for the LHC 
is unquestionable. 

Many efforts have been made in the last few 
years in developing efficient techniques for one- 
loop calculations, using both highly-refined tradi- 
tional methods and the new approaches of OPP 
and Generalized Unitarity. A significant achieve- 
ment of this work has been the calculation of 
complete NLO predictions for several interesting 



2^4 processes at LHC [2I3I4I5I6I7I8I9IT0J . 

In this contribution we describe the main fea- 
tures of the HELAC-NLD approach to NLO calcu- 
lations. We provide an overview of the compu- 
tational framework and briefly summarize some 
physics results obtained with this system. 

2. Structure of HELAC-NLO 

2.1. Virtual corrections 

The calculation of the virtual corrections is or- 
ganized in the framework of the OPP method. 

Any one-loop amplitude is characterized by a 
set of denominators given as products of loop 
propagators, Di{q) — (q +pt) 2 — mf, and a poly- 
nomial - in the loop momentum - function called 
numerator, Ni(q), such that it can be represented 
in the form 
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where both the loop momentum q and the numer- 
ator function arc considered to be evaluated in d 



dimensions. On the other hand, it is well known 
that any one-loop amplitude admits a decompo- 
sition in terms of a basis of d-dimensional scalar 
integrals. In the limit d — > 4, this decomposition 
can be cast into the form 
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where Box, Triangle, Bubble and Tadpole denote 
the known scalar loop functions and R is a ra- 
tional function in the invariants also known as 
rational term. 

The basic idea of OPP is to decompose the nu- 
merator of the integrand function, N(q), into a 
basis of monomials of the propagators appearing 
in the denominator: the coefficients of the ex- 
pansion are the a, &, c, d appearing in EqJ2] plus 
additional, ^-dependent, coefficients also known 
as spurious terms, whose analytical dependence 
on the loop momentum is known. It is therefore 
possible to evaluate the coefficients of the expan- 
sion by computing N(q) a sufficient number of 
times (equal the number of coefficients to be de- 
termined), by appropriate choices of q. One pos- 
sibility is to select those values of the loop mo- 
mentum that nullify certain sets of propagators 
- the so-called unitarity cuts - and solve the re- 
sulting linear system of equations in terms of the 
expansion coefficients. If the numerator N(q) is 
considered in d = 4 dimensions, it is necessary to 
decompose the rational terms, R — R± + i?2 : we 
call i?i the part of the rational terms that is de- 
rived through the reduction starting with a d = 4 
numerator function. i?2 is the part of the rational 
terms that originates from the explicit (d — 4) de- 
pendence of the numerator function evaluated in 
d dimensions. As has been proven in |12j i?2 can 
be fully determined by including new interaction- 
vertices with up to four legs, in full similarity with 
the usual UV counter-terms. 

On the one hand the method can be fully im- 
plemented numerically, resulting to a very fast re- 
duction of arbitrary multi-leg amplitudes. On the 
other hand, being inherently algebraic, provides 
exact results (with arbitrary precision) for the 
coefficients multiplying the one-loop scalar func- 



tions as well as for the rational terms. The OPP 
method thus provides a systematic framework to 
evaluate the cut-constructible part of the ampli- 
tude as well as the R\ term, given a numerator 
function as input (see |11I12| for details). 

In the approach we follow HELAC-1L00P [14], 
based on HELAC-PHEGAS [IS], automatically gen- 
erates all the input numerators for the OPP re- 
duction as well as all tree-order-like contributions 
including the new vertices needed for the full de- 
termination of the rational part (R2) |15I16| . as 
well as UV-renormalization counter-terms. For a 
detailed description of the algorithm we refer to 

El- 

To summarize, our approach for the evaluation 
of one-loop amplitudes relies upon three basic el- 
ements: 1) an implementation of the OPP reduc- 
tion; 2) a library of scalar loop integrals; 3) an 
algorithm to automatically evaluate numerator 
functions and R2 terms. These tasks are solved 
respectively by the codes CutTools, DneLDop and 
HELAC-1L00P J22I14J . 

We organize the calculation of virtual correc- 
tions using a re-weighting technique and adopt 
helicity and colour sampling methods in order to 
optimize the performance of our system [5] . 

The treatment of the QCD colour is based on 
the colour connection representation. Given a 

generic QCD amplitude, M.^"'?"' l '" ™ 9 , com- 
posed by n g gluons and n q quarks (antiquarks) 
with colour indices a and i (J) respectively, it 
is possible to get a uniform representation for 
quarks and gluons by contracting each gluon in- 
dex a with a corresponding £"■: 
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The resulting amplitude can be now decomposed 
in terms of colour-stripped partial amplitudes A a : 



M^;:::X = Y, 5 ^^^---^ k , 3k A„, (4) 



with k = n q + n g , and a denoting a permutation 
of the set {1, . . . , k}. Each A„ is interpreted as 
the contribution to the amplitude associated to 
one possible way of connecting the external par- 
ticles through colour flow, the connection being 
established by the indices of the <5's. The Feyn- 



man rules that allow the calculation of all A^s 
have been described in Refs. 23 24 13 . 

At the level of squared matrix element, the fol- 
lowing equivalence holds: 
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where the colour matrix, C a ^,, can be evaluated 
automatically by counting the number of cycles 
that are obtained by closing the colour flows of 
a, a'. Full colour summation is performed using 
the right-hand side of Eq. [5] 

While appealing, the colour-connection repre- 
sentation has a drawback in that the computa- 
tion becomes less efficient for processes with high 
parton multiplicities. Indeed, the number of par- 
tial amplitudes that have to be evaluated for a 
given phase-space point increases factorially with 
the number of the external coloured particles. A 
possible way out is given by the colour sampling 
approach where, following the left hand side of 
Eq. [5j one generates for each phase-space point 
a random colour assignment and computes the 
matrix element Ml 1 ' 1 , 2 ''"' 1 * accordingly: due to 
EqQJ only the subset of the A CT 's which have a 
colour flow compatible with the given assignment 
contribute. In fact, as shown in Table [TJ the num- 
ber of subamplitudes that have to be computed 
in this way is sensibly reduced, on the average, 
for increasing complexity of the colour structure. 

2.2. Real corrections 

The singularities arising from soft and collinear 
divergences in the initial and final states are 
treated within the framework of Catani-Seymour 
dipole subtraction 25i , in its formulation for mas- 
sive quarks |26| and arbitrary polarizations |27 . 
After combining virtual and real corrections, sin- 
gularities connected to collinear configurations in 
the final state as well as soft divergencies in the 
initial and final state automatically cancel for 
infrared-safe observables when the jet algorithm 
is applied. Singularities connected to collinear 
splittings in the initial state are removed via PDF 
factorization. 

Calculations are performed with the help of the 
code HELAC-DIPOLES [27] which implements au- 
tomatic construction of all relevant subprocesses 



and dipoles as well as phase-space integration of 
subtracted real radiation and of the correspond- 
ing integrated dipoles. The phase-space integra- 
tion is performed using the approaches of the 
multichannel generator PHEGAS [28, and of the re- 
cently published adaptive-multichannel generator 
KALEU [29] . 

A phase-space restriction is applied on the con- 
tribution of dipoles, as originally proposed in 
[30 31 , considering two values of the unphysical 
cutoff: a max = 1 and a max = 0.01. The a max in- 
dependence is explicitly checked in all our results 
as a consistency check of the calculation. 

3. Predictions for the LHC: tb + 2 jets 

The HELAC-NLO framework has been applied in 
the computation of NLO QCD corrections to sev- 
eral processes relevant for top quark and Higgs 
phenomenology [15 9. We report here, as an ex- 
ample, a brief summary of our results about inclu- 
sive production of it in association with two jets, 
which represents an important background for 
Higgs searches in the channel ttH(H —> bb) |32) . 

To be specific, we consider proton-proton colli- 
sions at the center of mass energy y/s = 14 TeV. 
Top quarks are considered on shell without any 
kinematical restriction applied. For the top quark 
mass we set m t = 172.6 GeV, whereas all the oher 
partons are considered massless. We consistently 
use the CTEQ6 set of parton distribution func- 
tions throughout the calculation, i.e. we take the 
CTEQ6L1 set with a one-loop running as at LO 
and CTEQ6M pdf 's with a two-loop running as 
for the evaluation of the NLO corrections. 

Our definition of jets is based on the kr algo- 
rithm [33134] using AR — 0.8, where AR denotes 
the separation in the azimuthal angle - pseudo- 
rapidity plane. Jets are obtained by recombining 
at most two partons and the recombination takes 
place only when both partons satisfy \yi\ < 5 (ap- 
proximate detector bounds). Also b quarks en- 
ter our jet definition, so the final state ttbb is in- 
cluded in the analysis. We further impose a min- 
imum AR =1.0 separation between jets and re- 
quire the latter to lie within the region defined by 
\yjet\ < 4.5. Finally, we require hard jets to have 
a minimum transverse momentum of 50 GeV. 



Table 1 

Comparison between the number of subamplitudes to be evaluated using full colour sum (n conn ) and 

using the Monte Carlo approach described in the text (average number, (n conn )Mc), f° r three reference 

processes. 



Process 



{nconn) MC 



gg -+ bbW+W- 


6 


1.74 


gg —> ttbb 


24 


3.04 


99 -> tigg 


120 


6.27 



Ratio 



3.5 

7.89 

19.1 



The numerical results of this study are summa- 
rized in Table [2] where we show a comparison of 
the total LO and NLO cross sections at the cen- 
tral scale [in = //j? = m t . As a cross check, we 
report the NLO results obtained for two different 
values of the dipoles cutoff parameter, a max = 1 
and a max = 0.01, showing that they are in mu- 
tual agreement within the integration error. 

The scale dependence of the corrections is illus- 
trated in Figure [TJ In the absence of any jet veto, 
QCD corrections determine a negative shift of 
approximately 11% compared to the LO results. 
A dramatic reduction of the scale uncertainty is 
observed when going from LO to NLO: varying 
the scale up and down by a factor 2 changes the 
LO cross section by +72% and -39% respectively, 
while in the NLO case one observes a variation of 
-13% and -12%. Applying a jet veto of 50 GeV 
results in an increasing negative shift operated by 
QCD corrections and leads to a K-factor of 0.64. 

Let us now analyse the changes in shape of dis- 
tributions determined by QCD corrections. A few 
differential cross sections for simple observables 
are considered. The most important one for us is 
the di-jet invariant mass since it is the observable 
entering Higgs boson studies. This is illustrated 
in Figure [21 where we focus the attention in the 
range up to 400 GeV due to our phenomenologi- 
cal motivations. We notice that the distribution 
shows tiny corrections up to at least 200 GeV. 

Much larger effects can be found in the other 
observables analysed, namely the transverse mo- 
mentum of the first and of the second hardest jet. 
Figures [3] and [4] show visible differences in shape 
in the region of high px, especially in the case of 
the first hardest jet. We conclude by pointing out 
that an accurate description of the tails of these 



Pt distributions cannot prescind from the inclu- 
sion of higher-order corrections. On the other 
hand, the behaviour at low px's is altered by soft 
collinear emissions whose effects can be simulated 
by parton showers. 




0.5 1.0 2.0 



Figure 1. Scale dependence of the total cross sec- 
tion for pp — »• ti + 2 jets + X at the LHC with 
Hr = I^f = £mo> where no = m t . The blue dot- 
ted curve corresponds to the LO, the red solid to 
the NLO result whereas the green dashed to the 
NLO results with a jet veto of 50 GeV. 



4. Outlook 

It is amazing how much progress has been re- 
alized over the last few years in the field of multi- 
leg NLO calculations. HELAC-NL0 is a complete 



Table 2 

Cross section at LO and NLO for ti+2 jets production at LHC (a/s = 14 TeV).The two NLO results 

reported refer to different values of the dipole phase-space cutoff a max . The scale choice is /i# = fip = m t . 



pp^tt + 2 jets + X 

no jet veto 

with jet veto (pr,x < 50 GeV) 



(Pb) 



,i (Pb) 



<CL=o.oi (Pb) 



120.17 ±0.08 
120.17 ±0.08 



106.95 ±0.17 
76.58 ±0.17 



106.56 ±0.31 
76.57 ±0.37 



io L 



> 

C3 



a. 



10" 



to 

-a 



10" 




100 



200 



300 



400 



mjj [GeV] 




100 200 300 400 

Pt.j [GeV] 



Figure 2. Distribution of the di-jet invariant mass 
for pp -> ti + 2 jets + X at the LHC. The red 
solid line refers to the NLO results while the blue 
dotted one to the LO one. 



Figure 3. Distributions of the transverse momen- 
tum for the first hardest jet for pp — > tt+2 jets+X 
at the LHC. The red solid line refers to the NLO 
results while the blue dotted one to the LO one. 



framework for NLO calculations, including vir- 
tual and real corrections. In its present form can 
be used for NLO QCD calculations for arbitrary 
processes including up to one-particle-irreducible 
8-point one-loop functions. It has been success- 
fully used for a number of processes of current 
interest, such as pp — > ttbb including the Higgs 
signal contributions and pp — > tt + 2 jets. In the 
immediate future, we plan a complete analysis of 
tt+jets production including the full description 
of top decays and non-resonant contributions, i.e. 
pp —> Ivl'v'bb. Plans of further development in- 
clude also the incorporation of all i?2 and UV 
counter-terms for full electroweak corrections. 
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